A&A manuscript no. 

(will be inserted by hand later ] 



Your thesaurus codes are: 
02.13.2;08.16.6;09.10.1 



ASTRONOMY 

AND 
ASTROPHYSICS 



On the physics of cold MHD winds from oblique rotators. 



on 

ON 
On 



3 

in 



in 
o 
l> 
o 

ON 
On 



Oh' 
i 1 

O 1 



X 



S.V.Bogovalov 

Moscow Engineering Physics Institute, Kashirskoje Shosse 31, Moscow, 115409, Russia, e-mail: bogoval@photon.mephi.ru 
Received date today; accepted date yesterday 



Abstract. I show that the self-consistent solution of the 
problem of MHD plasma flow in magnetosphere of an 
oblique rotator with an initially split-monopole magnetic 
field is reduced to the solution of the similar problem for 
the axisymmetric rotator. All properties of the MHD cold 
plasma flows from the axisymmetric rotators with the 
initially split-monopole magnetic field are valid for the 
oblique rotators as well. Rotational losses of the oblique 
rotator do not depend on the inclination angle and there 
is no temporal evolution of this angle. Self - consistent 
analytical and numerical solutions for the axisymmetric 
plasma flows obtained earlier show that the rotators can 
be divided on fast rotators ( (Jq/Uq > 1) and slow rota- 
tors ( <To/U§ < 1), where <7o is the ratio of the Poynting 
flux to the matter energy flux in the flow at the equator 
on the surface of the star, Uq = jqVq/c, vq and 70 being 
the initial velocity and the Lorentz-factor of the plasma. 
The self-consistent approximate analytical solution for the 
plasma flow from the oblique rotator is obtained under the 
condition cfq/Uq <C 1. Implications of these results for ra- 
dio pulsars are discussed. In particular, I argue that all 
radio pulsars are apparently the slow rotators ejecting the 
Poynting dominated relativistic wind. 

Key words: MHD - pulsars; jets and outflows 



1. Introduction 

An analysis of the relativistic plasma flow is necessary for 
understanding the processes taking place in radio pulsar's 
magnet osphe res, com pact g alactic objects and in AGN's 
(Arons |l996| ; Mirabel |l996] Pelletier et al. " 



199rj). In the 



present paper, as in previous one (Bogovalov 1997 ), we 
concentrate our attention totally on the problem of the 
relativistic plasma flow in the conditions typical for radio 
pulsars. 

In spite of systematic research in the field of the physics 
of radio pulsars in the last years, the structure of the 
magnetosphere and the mechanisms for the acceleration 
of plasma i n thes e objects to large extent remain vague 
(Lyubarsky 1995 ). One of the most important unsolved 



problem of the physics of radio pulsars is the problem of 
energetics of the ejected wind of relativistic plasma. For 
example, the kinetic energy of the plasma accelerated in 
the inner magnetosphere of the Crab pulsar is not suffi- 
cient to explain the energetics of the relativistic wind exit- 
ing the synchrotron nebula surrounding the pulsar (Arons 
1996b . 



The basis of the theory of the plasma production in 
radio pulsar magnetospheres was initiated by Sturrock 
(1971). This theory was developed in more detail for dif- 
ferent conditions on the stellar surface by Ruderman & 
Sutherland (1975) and by Arons ( [1981 ). Primary electrons 
are accelerated in the so called "electrostatic gaps" and 
produce dense relativistic plasma of secondary particles 
with Lorenz- factor 70 ~ 10 2 — 10 3 . This plasma screens 
the accelerating electric field and limits the potential drop 
by the value r yg ap mc 2 with -f gap ~ 10 7 . The total flux of 
the kinetic energy of all particles appears much less than 
the total rotational losses of the fast rotating radio pulsars 
due to this limitation. Almost all energy is carried out by 
the electromagnetic field. The wind with similar charac- 
teristics is form ed in the outer gap model (Cheng, Ho & 
Ruderman 1986| ). The wind of relativistic plasma from ra- 
dio pulsars can be characterized by the ratio of the Poynt- 
ing flux to the flux of the kinetic energy of the plasma a. 
The plasma is Poynting flux dominated when a > 1 and 
is kinetic energy dominated if a < 1. Electrostatic gaps 
give a > 1000 for the Crab pulsar. At the same time in- 
terpretation of observations of the Crab Nebula compels 
us to conclude that this ratio at large distances from the 
pulsar is very small. This conclusion is based on the as- 
sumption that the subsonic flow of the plasma terminated 
by the shock wave at the interaction of the wind with the 
interstellar medium can be considered as the flow of ideal 
plasma with the only dissipative process of synchrotron 
cooling. Under this assumption the observed expansion of 
the outer edge of the Crab Nebula, the synchrotron and 
TeV gamma-ray emission of the nebula produced via In- 
verse Compton Scattering of the relativistic particles on 
2.7-K MBR emission can be explain ed un der unique choice 
a = 3 ■ 10~ 3 (Atoyan & Aharonian |l996|) . 
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Recently Begelman ( |1998| ) has revisited the key as- 
sumption of the theory by Kennel & Coroniti (1984). He 
argued that even a wind with a ~ 1 in the pre-shock region 
can give the observable properties of the Crab Nebula if 
one takes into account dissipative processes in the nebula. 
It follows from the theory of relativistic MHD shocks that 
a wind with a ~ 1 in the pre-shock region creates the flow 
in the post-shock region with a > 1. But this flow must 
be strongly unstable. The instability provides an effective 
transformation of the magnetic field into the kinetic en- 
ergy of the plasma accompanied by acceleration of parti- 
cles in the nebula. It is important that in both (Kennel & 
Coroniti or Begelman) scenario it follows that there exists 
some unspecified mechanism for the transformation of the 
Poynting flux into the flux of the kinetic energy as the 
plasma travels from the star to infinity. This mechanism 
is apparently the basic mechanism for the plasma acceler- 
ation since it ensures the transformation of at least 50% 
( in the Begelman scenario) or 99.7 % ( in the Kennel & 
Coroniti model) of the rotational energy of the neutron 
star into the kinetic energy of the relativistic wind. 

One of the possible mechanisms of the acceleration is 
the magnetic acceleration of the plasma by the rotating 
magnetosphere. Unfortunately this mechanism appeared 
non effective for the axisymmetric rotators with typical 
pulsar's parameters. No effective acceleration was found 



There were a lot of attempts to solve this problem in 



neither in nearest zone ( Bogovalov 1997), nor in far zone 
of the rotator (Bogovalov & Tsinganos 1999). At the same 



time up to now there is no solution describing the MHD 
plasma flow from the oblique rotator. Usually it is be- 
lieved that the inclination of the magnetic moment of the 
star to the axis of rotation could be important for the 
plasma acceleration. In this paper we study this possibil- 
ity in the model of the oblique rotator with the initially 
split-monopole magnetic field on the surface of the star. 

2. Basic assumptions 

The axis of rotation of real radio pulsars is not directed 
along the magnetic moment. The solution of the prob- 
lem of the plasma flow in magnetosphere of this object 
is extremely complicated. To simplify the problem several 
models of radio pulsars were proposed. The sequence of 
these models is presented in Fig. |l|. Firstly Goldreich & Ju- 
lian (1969) proposed axisymmetrically rotating star with 
initially dipole magnetic field (step 1 in Fig. [I]). The ro- 
tational losses of the axisymmetrically rotating star eject- 
ing relativistic plasma are comparable with the rotational 
losses of the oblique dipole in vacuum. The energy of ro- 
tation is carried out near the surface of the axisymmetric 
rotator by the Poynting flux. This is why this model can 
be used for study of the process of the Poynting flux trans- 
formation into the kinetic energy of the plasma. 

However even with such simplification the structure 
of the axisymmetric flow of the plasma from the star with 
the initially dipole magnetic field appears too complicated. 



massless approximation (Michel 19731; Beskin, Gurevich 
& Istomin 1983; Lyubarskii 1990| ; Contopoulos, Kazanas 
& Fendt |1999[ ). To simplify the problem Michel Ql969| ) was 
the first who used the rotator with the prescribed split- 
monopole poloidal magnetic field for the investigation of 
the relativistic plasma flow in the magnetosphere of the 
axisymmetrical rotator (step 2). It follows from Fig. [jjthat 
there are no closed field lines in the split-monopole model. 
All lines go to infinity. This allows one to remarkably sim- 
plify the analysis of the flow. 



Michel (1969) considered the flow of the plasma in 
the prescribed split-monopole magnetic field. The affect 
of the moving plasma on the poloidal magnetic field was 
not taken into account. In the self-consistent solution the 
plasma and the electromagnetic field affect each other. 
The problem of the self-consistent plasma flow from the 
rotator with the initially split-monopole magnetic field 
was investigated in nonre lativi stic and relat ivistic limits 
in the paper s by Sakurai ( 1985| ), Bogo valov ( 1992 ) Bogo- 
valov ( |1997D , Bogovalov & Tsinganos (|l999| ). The phrase 
"initially split-monopole magnetic field" means that the 
normal component of the magnetic field on the surface 
of the star is constant but changes sign on the magnetic 
equator. In this paper we firstly consider the model of the 
oblique rotator with the split-monopole magnetic field on 
the surface of the star (step 3). This model allows one 
to investigate the plasma flow in conditions more typical 
for real radio pulsars than it occurs in the axisymmetric 
model and to connect the model with the split-monopole 
magnetic field with the real pulsars (step 4). 

Since the density of the relativistic plasma produced 
in the pulsar magnetosphere is high enough to screen the 
electric field parallel to the magnetic field, the magneto- 
hydrodynamical approximation is used in this paper to 
describe the flow of the plasma. The plasma is considered 



as cold to a first approximation (Lyubarsky 1995) 



3. Reduction of the oblique rotator problem to 
the axisymmetrical problem 

System of time dependent equations defining the temporal 
evolution of the flow of the relativistic plasma (Akhiezer 



et al 1975) 



(^T + ( vV )7v) =q-E+-j xH, 
at c 



<9H 

~c~di 



cur/E, 



curlll = — j + — , 

c cot 

divU = 0, 
div'Ei = 4:irq, 



(1) 
(2) 

(3) 

(4) 
(5) 
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Fig. 1. The sequence of models introduced to simplify the 
solution of the problem of the plasma outflows from radio 
pulsars. Step 1 is the transition from the oblique rotator 
with a dipolc magnetic field to the axisymmetric rotator. 
Step 2 is the transition from the axisymmetric rotator 
with the dipole magnetic field to the axisymmetric rotator 
with the split-monopole magnetic field. In this paper the 
transition from the axisymmetric split-monopole to the 
oblique split-monopole model is done (step 3) together 
with step 4 connecting this model with real pulsars. 



dn 
~dt 



+ div(vn) = 0. 



(6) 



together with frozen-in condition E|jvxH = gives in 
spherical coordinate system the equations of motion 
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and the induction equations 
8H r 1 
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cdt 
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(9) 



(10) 



(11) 



(12) 



The laws of conservation of the magnetic field and the 
matter flux, together with Coulomb law take a form 



1 d(r 2 H r 



+ 



1 d sin 9H g 



r 2 dr ' rsin# 

dn 1 d(r 2 nv r ) 1 
dt r 2 dr r siu 

1 d(r 2 E r ) 



+ 



1 dH^ 



d9 rsin6> dip 
d sin 9nvg 1 dnv, 



0. 



(13) 



d9 



dip 



= 0.(14) 



dr 



+ 



1 d sin 9Eg 



rsmf 



d9 



1 ggg 
rsin0 dip 



The frozen-in conditions are 



E r + -(v e H v - V v Hg) = 0, 



Eg + - {v v H r -v r H v ) = 0, 



E v + -(v r Hg - VgH r ) = 0. 



= Airq. (15) 

(16) 
(17) 
(18) 



In this system q is the induced space electric charge den- 
sity, 9 is the polar angle and ip is the azimuthal angle. 

Boundary conditions on the surface of the star with 
radius for the system of equations above are: 

1. The Lorenz- factor 70 is specified constant; 

2. The normal component of the density of the matter flux 
is specified constant and uniform; 

3. The normal component of the poloidal magnetic field 
Hq does not depend on coordinates and time in every point 
of the surface of the star and changes sign on the magnetic 
equator. 

4. The tangential component of the electric field is contin- 
uous on the star surface; 

It is assumed that the velocity of the plasma exceeds fast 
magnetosonic velocity (the flow is super sonic) at the infin- 
ity. This system of equations together with the boundary 
conditions describes the stationary axisymmetric flow as 
well. 

Let's assume that the self-consistent solution for the 
problem of the plasma flow from the axisymmetric rota- 
tor with an initially monopole-like poloidal magnetic field 
is specified. The phrase "initially monopole-like magnetic 
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Monopole-like magnetic field 



Splil-monopole magnetic field 



Fig. 2. The solution in the model with the split-monopole 
magnetic field can be obtained from the solution for the 
monopole-like magnetic field by change of sign of the mag- 
netic field in lower hemisphere. The same is valid for re- 
verse operation. There is a current sheet at the equatorial 
plane in the split-monopole model which is a contact dis- 
continuity in ideal MHD approximation 



field" means that the normal component of the poloidal 
magnetic field on the surface of the star is constant and 
does not change sign on the magnetic equator. Actually 
this field has no magnetic equator at all. Although such 
magnetic fields can not be created in reality, this rather 
artificial mathematical model is convenient to construct 
solutions for more realistic plasma flows. The monopole- 
like magnetic field differs from the split-monopole mag- 
netic field only by direction of the field li nes in one of 
hemispheres and was introduced by Michel ( 1969 ). In the 
model with the monopole-like magnetic field all field lines 
are out coming from the surface of the star as it is shown 
in Fig. ^. There is no current sheet in this flow in con- 
trast to the model with the split-monopole magnetic field 
which contains the current sheet in the equatorial plane. 
The current sheet is a contact discontinuity (Landau & 
Lifshitz 1963) in ideal MHD approximation. The magnetic 
field changes direction at the passage through the current 
sheet. It is obvious that the solution for the split-monopole 
can be obtained from the solution for the monopole-like 
magnetic filed by formal reverse sign of the magnetic fields 
in the lower hemisphere. Similar operation can be done in 
the case of the oblique rotator as well. 

The flow and the monopole-like poloidal magnetic 
field are perturbed by rotation in the self-consistent ax- 



isymmetric flow (Bogovalov & Tsinganos 1999). We as- 
sume that this perturbed self-consistent solution is known. 
Let's introduce the velocity V(r), the poloidal B p (r) and 
toroidal B v (r) magnetic fields and the density of the 
plasma N(r) in the self-consistent axisymmetric flow with 
the initially monopole-like magnetic field. The poloidal 
electric field E and the poloidal magnetic field are con- 
nected in the axisymmetric flow E = (r sin#Sl/c)B p x e v 
(Weber & Davis 1967), where e v is the unit vector di- 
rected in the azimuthal direction. The toroidal electric 
field is equal to 0. These variables depend only on coor- 
dinates. Dependence on time is absent for the stationary 
axisymmetric flow. This solution satisfies automatically 



all the equations of the system (|7T8) and corresponding 
boundary conditions. 

Now we show how the self-consistent solution for the 
plasma flow in the magnetosphere of the oblique rotator 
with the initially split - monopole magnetic field can be 
obtained from the known solution for the axisymmetric 
problem. Let's consider the following transformation of 
the axisymmetric solution. The velocity and the density 
are taken the same as for the axisymmetric rotator 



v(r,t) = V(r), n(r,t)=N(r). 



(19) 



The poloidal magnetic field H p (r, t) for the oblique ro- 
tator is obtained from the poloidal magnetic field of the 
axisymmetric rotator as follows 

H p (r,t)=?7(r,t)B p (r). (20) 

The same procedure gives us the toroidal magnetic field 

H v {r,t)=r,(r,t)B v {T), (21) 

where 7/(r, t) is unknown function to be specified. The 
poloidal electric field E is defined as in the axisymmet- 
ric case 

r sin 8Q 



En 



E r = r -^H 9 , 



(22) 



(23) 



and the toroidal electric field E v = 0. It is easy to show 
that this is indeed the solution of the problem for the 
oblique rotator and to specify the function rj. 

The frozen-in conditions ( 16-18]) are fulfilled automat- 
ically for the solution (|l9|-^3|). At the stationary rotation 
the dependence of rj on ip and t can be presented by one 
variable £ = (p — Qt in the spherical system of coordinates. 
Thus, rj depends on three variables r, 6 and £. Then the 
induction equation ( |Io|) is reduced to the equation 



d E„ 



{111, 



-) = o 



(24) 



<9£ r sin ( 

which is fulfilled due to condition (p2|). The same is valid 
for the induction equation (^) which is reduced to the 
equation 



3_ E, 



dt; r sin 6 



-) = o 



(25) 



which is satisfied due to condition (|23|). Third induction 
equation (|l2|) after substitution of conditions ( |22] , p3| ) is re- 
duced to the flux conservation condition (|l3|) which takes 
a form 



(B • Vt?) = 0. 



(26) 



since the total magnetic field B of the axisymmetric flow 
satisfies the equation V • B = 0. It implies that function 
i] is constant on the field line of the magnetic filed of the 
axisymmetric flow. 
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Fig. 3. The initial position of the star and the geometry- 
used in equation 



Fig. H shows the initial position of the star with the 
magnetic moment M inclined to the axis of rotation at 
the angle x at i = 0. It follows from equation (|2(]) and the 
boundary conditions that the function rj is equal to 1 in 
the points on the surface of the star where the field line 
go out from the surface of the star and is equal to —1 in 
the points where the filed lines come into the surface. It 
is convenient to introduce function 



D(x) = 



1, if x > 1 
-1, if x < 1. 



(27) 



The sign of the function r\ on the surface of the star varies 
with the sign of the product (e • e^/), where e is the unit 
vector directed to the point on the surface of the star, 
is the unit vector directed along the magnetic moment. 
This product can be presented as 

(e • ej,f) = sin x sin 9 sin ip + cos \ cos 9. (28) 

Then, on the surface of the star, the function r\ is 

rj(9, ip — fit) = D (sin x sin 9 s'm((p — fit) + cos x cos 9). (29) 

Actually this is the boundary condition to equation (p(f). 
The equation can be presented in the form 



dr) dr/ drj 

Br ~x~ + Be ^a + B f — ■ ax 
or ro9 r sin 9d(p 



0. 



Equations of characteristics for this equation are 
dr r sin 9dip 



B r 



I' 



and 

dr rd9 

Therefore the general solution is 
r r B H dr 



r,(r,e,<p,t) = f(9 



rB r 



(30) 



(31) 



(32) 



^%r\ (33) 
• sin 9B r 



where the integrals over r are taken along the field line of 
the poloidal magnetic field of the axisymmetric solution, 
/ is arbitrary function. The solution satisfying boundary 
condition (^9|) has a form 

r)(r,9,(p,t) = D(sm(x)am(9- J 

x M<P ~ Ik 7§^k - ™) + 
+ ™s(9-f R ^)cosx 



(34) 



It follows from this solution that rj 2 = 1 and rj changes 
sign when the magnetic field changes direction. It is easy 
to show now that equations of motion (fij|9|) are also sat- 
isfied for the solution jis|-p3|) at the function r\ defined by 
([34j). Notice that on the left hand side of these equations 
there is no function r\. In the right hand side of the equa- 
tions of motion function r\ comes in combination Aii] d ^ k , 
where A4 and B^ are arbitrary components of fields of the 
axisymmetric solution, xi spatial or time coordinates in 
4-space. This relationship can be presented as 



dxi 



dx. 



A i B k 1 -^ = A i ^. (35) 
2 dxi dxi ' 



Therefore, the function 77 disappears in the equations of 
motion. Here we ignore the difference in the dynamics of 
the current sheet and the surrounding plasma assuming 
that the current sheet is the mathematical discontinuity 
as usual in ideal MHD. This assumption can be violated 
for the oblique rotators at large distance from the star. But 
at large distances the dynamics of the current sheet can be 
considered particularly in WKB approximation (Coroniti 
|1990| ; Michel |l994|) . 

Thus, we obtain the self-consistent solution for the 
oblique rotator from the known self-consistent solution 
for the axisymmetric rotator. The sketch demonstrating 
the structure of the cold wind from the oblique rotator 
is presented in Fig. ||. The structure of the plasma flow 
is symmetric in relation to the equator. The form of the 
poloidal field lines is the same as for the axisymmetric ro- 
tator. In general there is collimation of the plasma flow 
to the axis of rotation, although the effect of the collima- 
tion depends on the parameters of the problem (Bogovalov 
& Tsinganos 1999] ). In the axisymmetric flow the current 
sheet dividing the magnetic fluxes of opposite directions 
is located on the equator. In the wind from the oblique 
rotator the current sheet takes a form of a wave. In the 
poloidal plane the poloidal magnetic field lines change di- 
rection on the current sheet. At first glance it seems that 
this behavior contradicts to the magnetic field freezing. 
The bottom panel of Fig.|] shows the structure of the field 
lines in the equatorial plane. It is seen that there is no 
contradiction with the magnetic flux freezing since the to- 
tal magnetic field depends on the azimuthal angle ip. The 
velocity and the density of the plasma do not depend on x 
and ip and are the same as for the axisymmetric rotator. 
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is denned by torques on the surface of the star (Michel & 
Goldwire p970|) 
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Fig. 4. Top panel shows the structure of field lines and the 
current sheet (thick wave-like line) in the poloidal plane. 
The lower panel shows the same in the equatorial plane. 
Arrows show the direction of the magnetic field lines. The 
direction of the field lines changes on the current sheet. 



Only the magnetic and the electric fields are modulated 
with the period of rotation. The sign of these fields is mod- 
ulated in the sector limited by the field lines which have 
roots on the magnetic equator. The squares of these fields 
are not modulated. The wave is simply the wave of the 
contact discontinuity (Landau & Lifshitz 1963 ) which is 
convected with the velocity of the plasma. 

The rotational losses of the oblique rotator can be cal- 
culated as the total Poynting flux through a sphere sur- 
rounding the star. The Poynting flux is defined by the 



formula c ^^ 1 . This formula is biquadratic on the com- 
ponents of the electromagnetic field. Since for the oblique 
rotator these components differ from the component for 
the axisymmetric rotator only on function 77, the Poynting 
flux appears independent on the function r\ and therefore 
the rotational losses are independent on the angle \. 

There is also no alignment of the magnetic moment 
along the axis of rotation as it happens in vacuum ap- 
proximation ( Ostriker & Gunn 1969] , Pacini 1965, Michel 



& Goldwire PL970 



Davis & Goldstein 197C) or disalign- 
ment as i t hap pens in the solution by Beskin, Gurevich & 
Istomin ( |l983| ). The evolution of the inclination angle \ 



Ly = 



L, 



[T r g COS ip ■ 



Trip cos 9 sin <p] d 2 S, 



[T r g cos <p — T rip cos 9 sin ip] d 2 S, 



and 

L z = / T rv sin9 d 2 S, , 



(36) 



(37) 



(38) 



where T r g and T rtp are the components of the energy- 
momentum tensor T r g and integration in ( p7| ) is performed 
over the surface of the star. The energy-momentum tensor 
is also biquadratic on the components of the electromag- 
netic field (Landau & Lifshitz 1975). Therefore it does 
not depend on the function rj and the angles x an d <P- The 
torques for the oblique rotator are the same as for the 
axisymmetric rotator. Integration in eq. (37) over the sur- 
face of the star gives L x = L y — 0. There are no torques 
which can change the inclination angle of the oblique rota- 
tor with the initially split-monopole magnetic field. This 
is not totally unexpected conclusion. The torques on the 
star rotating in vacuum differs from the torque on the 
start ejecting ideal plasma. The torque on the star rotat- 
ing in vacuum always align the rotational and magnetic 
axes (Sopcr 1972). The temporal evolution of the inclina- 
tion angle \ of the star ejecting plasma depends on the 
distribution of the magnetic field on the surface of the 
star. It was shown by perturbation theory in first order 
approximation on f2 that magnetic and rotation axes tend 
to align when the magnetic flux is more concentrated to 
the magnetic poles and evolve to x = 90° when there is 
less flux density at the poles than at the equator (Mestel & 
Selley 1970). Thus, the result obtained here for arbitrary 
angular velocity of the rotator with the split-monopole 
magnetic field is due to combination of two conditions: by 
the ejection of the plasma simultaneously with the uni- 
form magnetic flux distribution on the surface of the star. 

4. The flow of the plasma at < 1. 

The solution of the problem of the cold plasma flow from 
the rotator with the initially split-monopole magnetic field 
is defined by the following parameters and variables: H , 
mno, vq, fl, c, r, 9, <p, R*,x an d t. Lower index "0" denotes 
the values on the surface of the star. Since the dependence 
on x, V> an d t is known, we can consider only the depen- 
dence of the solution for the axisymmetric flow on other 8 
parameters. There are 3 parameters in this set with inde- 
pendent dimension. Therefore, according to the theory of 
dimensions (Barenblatt 197SQ the solution in dimension- 
less variables depends only on 5 parameters. Let us write 
this dependence as 



A = f(R star /a, al, a2, r/a, 



(39) 
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Here A is the vector of the state of the plasma includ- 
ing density, velocity and magnetic field, f is the unknown 
vector-function and a is some parameter with dimension 
of length. It is reasonable to consider limiting case when 
the dimensionless radius of the star goes to zero. It allows 
one to decrease the number of variables in the solution 
on one. Now the solution crucially depends only on two 
parameters and two variables in the dimensionless presen- 
tation. 

Let's measure the velocity of the plasma in units of the 
initial velocity vq and all geometrical variables measure 
in the initial radius of the fast mode surface (FMS) r/. 
The FMS is the surface where the poloidal velocity of the 
plasma equals the local velocity of the fast mode MHD 
perturbations. In the comoving coordinate system where 
the electric field is equal to zero, this condition takes a 
form 



H 



■litmn c 2 + H" 2 ' 



(40) 



where n and H are the density of the plasma and the 
magnetic filed in the comoving system. Taking into ac- 
count that H 2 — E 2 is invariant (Landau & Lifshitz 1975| ) , 
77 can be calculated at f2 = as follows 



TJ2 E>2 



(41) 



f — ■ 

1 AirmcnoVoUo 

In these units the radius- vector is X = r/r/. The electric 
and magnetic fields can be measured in the units of the 
initial poloidal magnetic field on the FMS. The density 
of the plasma can also be measured in units of the initial 
density on the FMS. 

The parameters al and a2 can be chosen rather arbi- 
trary. In particular it is convenient to choose al = a = 
(r/fi/uo7o) and a2 = e = (r/O/c). Let's consider the de- 
pendence of the solution on these parameters. 

Nonrclativistic limit corresponds e -t and vq <C c. 
The solution depends only on one parameter in this limit. 
The axisymmetrical nonrelativistic plasma flow was inves- 
tigated numerically by Bogovalov & Tsinganos (1999) for 
the parameter a achieving the value up to 4.5. It was found 
that there is clear division of the rotators on the fast and 
slow in dependence on the value of a. The slow rotators 
correspond to a < 1 and the fast rotators correspond to 
a > 1. There is strong physical difference between them. 
The plasma in outflows from the slow rotators is not accel- 
erated considerably and is not collimated in the subsonic 
region. In the case of the fast rotators the plasma is effec- 
tively accelerated and collimated to the axis of rotation 
already in the subsonic region of the flow. 

The flow of the relativistic plasma depends crucially on 
two parameters. But in this case it is also possible to divide 
the rotators on the slow and the fast ones depending on 
value of the parameter a. The solution can be presented as 
a vector function A with components A\ which are density 
of plasma, components of velocity, magnetic field and etc. 



Let's consider the expansion of this function on powers of 
a in the point a — at fixed e ^ 0. Physically it means 
that we consider the relativistic limit of the problem at 
70 — > 00. The expansion can be presented as 

A, = // 0) (6, X)(l + af W (e, X) + (afi ^ (e, X) + ....). (42) 

The first term /> (e, X) in this expansion has very sim- 
ple presentation at arbitrary parameter e ^ 0. It was 
shown that at the limit under consideration the function 
fi°\e,X) is following in the ordinary (dimension) vari- 
ables (Bogovalov, 1997) 



U r = U , U V =Q, Ug=0, n{r) 



no(-) ■ 



and 



, ril sin 6 , 



)H p , E 



Ho(^) 2 , if*< tt/2 
-Ho(^) 2 , i£6> tt/2 



H e = 0. 



(43) 



(44) 



(45) 



This is exactly the solution obtained by Michel (1973a) for 



the massless plasma. It is easy to show that the full system 
of the equations of motion is satisfied for this solution 
except the frozen- in condition ([l8|). The residual in this 
condition on the FMS is 

, («A) 2 



6 v, 



(1 + vq/c)' 



(46) 



The residual shows that the corrections to the solution are 
indeed proportional to powers of a and can be neglected in 
the region limited by the FMS provided that a C 1. This 
residual also shows that expansion (|42) contains terms 
(a/e) k , where k is integer number. Therefore the first term 
in expansion ( |42"| ) gives incorrect result in the limit a — > 
0, e — > provided that a/e = const which corresponds 
to the slow rotation of the star at arbitrary velocity vq. 
Mathematically it means that the point (a = 0, e = 0) is 
a particular point of the expansion. The terms (e/a) k can 
be eliminated from expansion (|42|) if to define the function 



gi(a/e,X) 



lim 

a/e = const 



Ai(a,e,X) 



(0) 



(e,X) 



(47) 



Then the expansion can be presented as 

A l = // 0) (e, X)(g(a/e, X) + £ a k g[ k) (e, X)) 

fc=i 



(48) 



The expansion X}fc=i ak 9^ ( £ j X) already does not con- 
tain terms (a/e) k and function Fi = f[(e,X)(g(a/e,X) 
gives correct first term at a — > for arbitrary parameter 
e including point e = 0. This function was calculated an- 
alytically by Bogovalov (1992). The corrected solution of 
the problem in the limit a — > at arbitrary e coincides 
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with the solution (43-45|) with exception of the toroidal 
magnetic field. It is replaced by 



(49) 



According to this solution the plasma moves radially 
in the split-monopole poloidal magnetic field with con- 
stant poloidal velocity and without any toroidal motion. 
At large distances the collimation of the plasma to the axis 
of rotation and corresponding acceleration of the plasma 
are also very weak (Bogovalov & Tsinganos 1999). It is 
reasonable to call all rotators with a < 1 as the slow ro- 
tators and the rotators with a > 1 as the fast rotators 
independently on the value of e. 

Let's consider the physical sense of the parameters a 
and e. It is clear that e is the ratio of the initial radius of 
the fast mode surface to the light cylinder. On the other 
hand 



e = oo = 



H 2 



47mom,c 2 7o vq 



(50) 



For the slow rotators <ro is the ratio of the Poynting flux 
to the matter energy flux (kinetic energy flux for the rela- 
tivistic plasma) at the equator and coincides with the well 



known magnetization parameter (Arons 1996) for the slow 
rotators. It is easy to show that 



op 

U 2 ' 
L o 



(51) 



Therefore the slow rotators correspond to the condition 
p « 1. 

n 

Diagram in Fig. H shows different regimes of the flow 
from rotators with the initially split-monopole magnetic 
field in coordinates <jq an d Uq. The line <tq = 1 divides the 
rotators with the Poynting flux dominated flow near the 
surface of the star (above this line) from the rotators with 
the matter energy dominated flow (below this line). The 
slow rotators are located below the line ero = Uq while the 
fast rotators are located above this line. In wide range of 
the parameters the regimes of the axisymmetric flows have 
been already investigated. Now it is evident that the ob- 
tained solutions are valid for the oblique rotators also. For 
the slow rotators with the strongly matter energy domi- 
nated flow (ero — > 0) an approximate analytical solution 
with accuracy to the terms of the order a 2 was obtained 



by Bogovalov (1992). At the increase of oo at constant 
Uo, rotators with the nonrelativistic flow firstly cross the 
line ero = Uq and becomes fast rotators. Numerical solu- 
tions for the nonrelativistic plasma flows (uq — > 0) from 
fast rotators were obtained up to a — 4.5 (Bogovalov & 
Tsinganos [19991 ). It was found that due to the strong col- 
limation of the plasma to the axis of rotation, the flow be- 
comes subsonic in some region along the axis of rotation 
and therefore should be absolutely unstable well before 
t he ro tator reaches th e line <jq = 1 (see also Lyubarskii 
(|1992D and Begelman (|l99S|) ). 
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Fig. 5. Regimes of the cold plasma flows. The locations 
of 6 pulsars emitting observable gamma-rays above 100 
MeV: Crab -1, PSR 1508-59 - 2, Vela - 3, PSR 1706-44 
- 4, PSR 1951+32 - 5, Geminga -6 and PSR 1055-52 - 7 
are shown by dashed lines for all possible values of Uo- 
All the pulsars are located in the region of slow rotators 
ejecting the Poynting flux dominated wind provided that 
Uq > 200. 



The rotators with the initially relativistic plasma flow 
firstly cross the line ctq = 1 at the increase of ao at 
constant Uq and became the slow rotators ejecting the 
Poynting flux dominated wind. At the increase of a they 
also cross the line a = Uq and become the fast rotators. 
The relativistic stationary plasma flow from the fast ro- 
tators was investigated recently by Beskin, Kuznetsova 
& Rafikov ( 1998j ). Surprisingly, the authors found that 
the perturbation of the poloidal magnetic field go to 
zero as ctq/Uq — ► and acceleration of the plasma re- 
mains very non effective like in the model by Michel 
( 1969[ ) with the prescribed monopole-like poloidal mag- 
netic field. The Lorentz-factor of the plasma reaches the 
value 700 ~ (er 7o) 1/3 . 

The plasma flow from the slow oblique rotator becomes 
especially simple. In particular, the function rj for the slow 
rotators (a <C 1) is 



rj(r, 6*, ip, t) = D(sinxsm0sin(<^ — Cl(t — r/vo))+ 

+ cos 8 cos x) . (52) 

Rotational losses of the slow rotators with the ini- 
tially split-monopole magnetic field are defined by the to- 
tal Poynting flux through the surface of the star 



2 H 2 Rjn 2 

3 Vq 



(53) 
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and they do not depend on the angle x- It is useful also to 
express the rotational losses through the total magnetic 
flux t/j of the open field lines of one direction (for example 
outcoming from the star) 



E r „t — 



1 V 2 Q 2 



(54) 



6tt 2 vq 

where ip = 2irH Rl for the split-monopole magnetic field. 

5. Implications to radio pulsars. 

How the model considered above is connected with the real 
radio pulsars? It is reasonable to assume that the rotator 
with the split-monopole magnetic field is equivalent to real 
radio pulsar if it provides the plasma flow similar to the 
plasma flow from the real pulsar at distances larger than 
the light cylinder. Therefore the equivalent split-monopole 
rotator should have the matter flux, the poloidal magnetic 
field flux from one of the poles, 70, vq, inclination angle x 
and the same as the real pulsar has. 

The magnetic flux from one of the poles of the split- 
monopole rotator is ip — 2ttR^Hq. The magnetic flux of 
the open field lines from the radio pulsar with dipole 
magnetic field is defined by the radius of the polar cap 

R p = R*\j~^p- (Goldreich & Julian |l969| ) . Therefore the 
flux of the open field lines from one of the polar caps of 
the radio pulsar is estimated as 

ijj = ttHoRK — ). (55) 
c 

The ratio of the poloidal magnetic field to the matter flux 
density H p /nv p is constant on the field line. Therefore 
this ratio for the split-monopole rotator and for the radio 
pulsar should be the same. These conditions allow one to 
define all parameters. The split-monopole rotator equiva- 
lent to the pulsar with the magnetic field on the polar cap 
Hq, angular velocity Q and ejecting plasma with initial 
density no and 70 has the magnetization parameter 

Hi RSI o 

<?r = ~ ^ ( — )*. (56) 



87rnomc 2 7o ' 

Below we consider the relativistic plasma with vq = c. The 
rotational losses of this object are 



1 H*Rltf 



6 c3 < 57 > 
These rotational losses equal the losses of the dipole ro- 
tating in vacuum at the angle of inclination x — 30°. 

The initial density of the plasma no is defined by 
the processes of multiplication of plasma in the electro- 
magnetic cascades initiated by primary particles with the 
Lorentz-factor 7 gop ~ 2 • 10 7 in the magnetosphere of pul- 
sar. In the inner gap theories (Ruderman & Sutherland 
1975; Arons |1981| ) the primary beam has Goldreich- Julian 



density 



nH 

2irec 



(58) 



We neglect here the dependence of uqj on the inclination 
angle. The primary particles produce A secondary elec- 
trons and positrons. Therefore the initial density of the 
plasma is no = \uqj. Calculations of cascades performed 



by Daugherty & Harding (1982) and by Gurevich & Is- 
tomin ( |1985| ) show that A ~ 10 4 . Outer gap model (Cheng 
& Ruderman |1986 ; Romani 1996|) g ives similar values of 
no. After substitution of (|5"8|) in (|56|), a r takes a form 



eHn 



^ RS ^ 3 



(59) 



The estimates of 70 are less definite. The flow of the 
plasma formed in the electromagnetic cascade consists of 
two components. One of them is the beam of the primary 
particles which can loss remarkable amount of energy. An- 
other component is the plasma of the secondary particles. 
It is reasonable to take the average Lorentz-factor of the 
secondary particles as 70 because it is this plasma which 
screens the electric field and ensures the frozen-in condi- 
tion. The average Lorentz-factor of the secondary particles 
strongly depends on the model of the gap. This parame- 
ter lies in the range 20 < 70 < 7 ga p/A. To be definite we 
assume here that 



70 = 0.1 



A 



200. 



(60) 



This value well agrees with the results of cascade simu- 
lations performed by Daugherty & Harding ( |l982 ). For 
other 70 the magnetization parameter is scaled as o> = 
0200 (2OO/70). The location of all radio pulsar observed 
by EGRET in gamma-rays above 100 MeV (Crab, Vela, 
Geminga, PSR 1706-44, PSR 1509-58, PSR 1951+32 and 
PSR 1055-55; Thompson 1998 ) in the parameter space ay, 
Uq for all possible Uo is shown in Fig. ^ by dashed lines. For 
the estimates of the magnetic field Hq equation ( |57j ) was 
used. It was assumed also that the momentum of inertia 
of all pulsars is 10 45 g • cm 2 . It is seen that all the pulsars 
appears the slow rotators provided that Uq > 200. This 
implies that they do not accelerate plasma and do not col- 
limate it considerably to the axis of rotation beyond the 
light cylinder if the flow of the plasma is dissipativeless 



(Bogovalov 1998; Bogovalov & Tsinganos 1999) 



6. Conclusion 

This paper clearly demonstrates that the oblique rotators 
do not differ strongly from the axisymmetric rotators. In 
the particular case considered in this paper the dynam- 
ics of the plasma ejected by the oblique rotator is exactly 
the same as the plasma dynamics of the wind from the 
axisymmetric rotator. It allows us to apply the all results 
obtained for the cold axisymmetric plasma flows to the 
oblique rotators. But these results show that for parame- 
ters typical for radio pulsars the acceleration of the rela- 
tivistic plasma is not effective. Radio pulsars apparently 
are basically slow rotators. The obliqueness of the rotators 
does not change this conclusion. 
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Our results indicate that only violation of the ideal 
MHD can provide effective acceleration of the winds from 
radio pulsars since neither magneto-centrifugal mecha- 
nism or the acceleration by the gradient of the toroidal 
magnetic field at the collimation of the plasma to the axis 
of rotation are effective for the pulsars. The role of the 
dissipative processes in the relativistic winds is under dis- 



cussion for many years (Coroniti 199C, Michel 1994, Mestel 



& Shibata 1994). Recently Melatos & Melrose (1996) tried 



to argue that the violation of the frozen-in condition takes 
place in the wind from an oblique rotator and it can pro- 
vide the acceleration of the wind. Our self-consistent solu- 
tion of the problem for the particular case of the oblique 
rotator with the initially split-monopole magnetic field 
shows no evidence of the violation of the ideal MHD in 
the wind. 



Mestel & Shibata ( 1994 ) argued that the frozen-in con- 
dition can be violated right after the light cylinder due to 
the formation of singularities in the stationary axisym- 
metric flow. This assumption is not confirmed by last re- 



sults by Bogovalov (1997) and by Contopoulos, Kazanas 
& Fendt dl999D . 

It follows from all results obtained recently in the 
physics of relativistic plasma outflows that only dissipa- 
tive processes can ensure plasma acceleration under pul- 
sar conditions. One of these mechanisms was proposed by 
Coroniti Q199CD and Michel fll994Q . They argued that due 
to the geometrical reasons the dissipative processes can 
ensure the transformation of the Poynting flux into the 
kinetic energy of the plasma provided that the inclination 
angle is high enough (~ 80°). The Poynting flux is trans- 
formed into the thermal energy of the plasma which is 
than accelerated by the thermal pressure gradient. 

Another possibility is the Poynting flux transformation 
in the volume of the flow. This can happen if the Poynting 
flux dominated plasma flow is unstable in relation to MHD 
perturbations. The anomalously low conductivity due to 
the MHD turbulence would provide the violation of the 
ideality in the flow with corresponding transformation of 
the Poynting flux into the energy of the plasma. The ad- 
vantage of this approach is that it can provide effective 
acceleration even at the axisymmetric rotation. But the 
question about stability of the Poynting flux dominated 
plasma flow is still open and should be investigated in 
details separatly. 
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